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We consider the problem of estimating the phase of squeezed vacuum states within a Bayesian 
framework. We derive bounds on the average Holevo variance for an arbitrary number A'^ of uncorre- 
lated copies. We find that it scales with the mean photon number, n, as dictated by the Heisenberg 
limit, i.e., as n~^, only for N > 4. For N < 4 this fundamental scaling breaks down and it be- 



comes n 



-N/2 



Thus, a single squeezed vacuum state performs worse than a single coherent state 
with the same energy. We find the optimal splitting of a fixed given energy among various copies. 
We also compute the variance for repeated individual measurements (without classical communica- 
tion or adaptivity) and find that the standard Heisenberg-Iimited scaling n~^ is recovered for large 
samples. 
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I. INTRODUCTION 

Squeezed states can improve the sensitivity of laser in- 
tcrferometry with a precision such as to beat the so-called 
shot-noise limit [1] and, therefore, have been considered 
as useful states in a wide variety of applications. They 
offer an enhanced resolution-energy tradeoff as compared 
to coherent state intcrferometry [2-4] . Some recent appli- 
cations include precision measurements of distances [5] , 
detection of small displacements in optical images [6], 
or optical imaging [7] with multimode light. One of the 
most promising fields of application is the detection of 
gravitational waves. This idea has been recurrently dis- 
cussed in the last two decades [4] and it is finally being 
included in the latest experimental proposals [8-10]. Ex- 
tensions to non-optical systems such as squeezed atomic 
states have also been considered [11-13]. 

The main advantage of squeezed states is that they can 
have optical phase variance below the standard quan- 
tum limit. In this context the most relevant ones are 
the squeezed vacuum states (SVS) [14]. It is generally 
claimed that these states have a phase variance scaling 
as n~^, where n is their average photon number. This 
result is hinted at by a somewhat heuristic argument that 
uses the Heisenberg uncertainty relation, and so, this 
n~^ scaling is known as the Heisenberg Limit (HL). A 
rigorous proof was derived by Holevo [15, 16], and can 
also be obtained from the Cramer-Rao bound [17] and 
the Braunstcin-Caves information inequality [18]. This 
work served as the foundation for recent developments in 
entanglement-enhanced metrology [19, 20]. The bounds 
on the variance derived in [14-20], however, are generally 
tight only when a large number of independent trials are 
repeated, whereas for small number of trials the attain- 
ability is not guaranteed. 

Our aim here is to determine the attainable precision 
in phase estimation for an arbitrary number N, not nec- 
essarily large, of uncorrelated, identically prepared SVS. 
To the best of our knowledge, no study has yet addressed 
this issue. Our analysis may be relevant in situations of 



quantum-limited communication in which one has access 
to few states (see [21]) or in gravitational wave detec- 
tion, where the available exposure time is limited and 
a large sample cannot be assumed. To tackle this situa- 
tion we adopt a Bayesian approach, in which an averaged 
cost function is minimized over all possible estimation 
strategies. For the problem at hand, the Holevo phase 
variance [15] is particularly well suited. This approach 
will enable us to perform a complete analytical computa- 
tion and derive closed expressions for the optimal aver- 
aged phase variance. We pay particular attention to the 
asymptotic values of this variance for SVS with large and 
small mean photon number. We also study the large N 
regime and compare the values with those for repeated 
individual measurements. 

We obtain the rather surprising result that the scal- 
ing n~^ of the Heisenberg limited variance of A'' SVS 
can only be attained if N is greater than four, while for 
1 < iV < 4 the scaling is given by n~^/^ (see Fig. 1 be- 
low) . This shows that a single SVS cannot do better than 
a single coherent state with the same energy. For large N, 
we, of course, recover the HL expression. Our analysis 
also shows that for a given amount of energy E that can 
be split among N identical SVS states of mean photon 
number n, i.e. E = nN (in units of hcu), the minimum 
variance is not attained with a single state, A^ = 1, of 
maximal mean photon number, as the n~^ scaling would 
suggest. For a large amount of available energy, we find 
that the optimal choice consists of splitting the energy 
among A^ = 8 identical copies (see Fig. 2 in Sec. HI). 

In this work we consider the most general measure- 
ments for optimizing the resolution and, therefore, our 
results can be regarded as the maximum precision limits 
to phase estimation with SVS allowed by quantum me- 
chanics. In this sense, our results are also relevant for 
metrology. Although there exist some quantum metrol- 
ogy precision bounds [20] that can surpass the HL, they 
require the use of nonlinear Hamiltonians, which arc in 
general very difficult to implement. 

The paper is organized as follows. In the next section 
we present the basics of our approach by considering the 



simpler case of coherent states. This serves as a warm- 
up exercise and as an illustration of the main techniques 
used throughout this paper. In Section III wc move to 
the core of our work. We derive the optimal measurement 
for an arbitrary number of identical copies of SVS and wc 
obtain the bound on the phase variance. In Section IV 
we provide results for non-adaptive individual measure- 
ments. In Section V wc discuss the implications of our 
results and draw our conclusions. The paper ends with 
two technical appendices. 



II. COHERENT STATES 

In this section we introduce the basic elements of 
our approach and consider the simpler case of coher- 
ent states as an example. These states are given by 
\tp) = £)(a)|0), where D{a) is the displacement oper- 
ator D{a.) = exp[aa'^ — a* a] and a [a^) are the stan- 
dard photon anihilation (creation) operators of a generic 
mode. The state of N identical copies of \4') is simply 
1$) ^ l^/))®^ or, equivalently, p = |$)($|. The phase 
shifted state \'4>0) is the result of acting on |-0) with 
the unitary operator U{0) = exp[i0a^a] (e.g., the time 
evolution operator for the free electromagnetic mode: 
= hut), i.e., IV'e) = U{0)\^p). The set of N identi- 
cal shifted states are likewise written as |$g) = \ipe)'^^ 
OT pe = |$e)($e|- 

Since our purpose is to estimate 9 with highest preci- 
sion, we allow ourselves to perform the most general mea- 
surements on \^e)- These, so-called generalized measure- 
ments, are described by a Positive Operator Valued Mea- 
sure (POVM); that is, a set Al = {Ox} oi positive oper- 
ators, O^ > 0, that add up to the identity, X)v '^x ~ ^■ 
Also, for each outcome of the measurement we need to 
give a guess of the value oi9, in technical words, we have 
to choose an appropriate estimator 9 : x ^'^ ^^x- ^^ ^^^ 
Bayesian approach this is done by providing a cost func- 
tion (figure of merit) for the estimation protocol. The op- 
timal choice of both measurement and estimator is taken 
to be that that minimizes (maximizes) the averaged cost 
function (figure of merit). We take the cost function 



V = Ue'^"-'^) 



1, 



(1) 



which is the natural variance for cyclic variables, as 
pointed out by Holevo [15]. Note that V vanishes for 
perfect estimation, and goes to infinity for a flat distri- 
bution of 9 (random guessing). Notice also that the min- 
imization of the variance V does not guarantee that 9 is 
close to 9 unless some sort of unbiasedness condition is 
imposed. For the distributions in the circle considered 
here the following conditions are sufficient: 
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We find it also convenient to introduce the figure of merit 
I , which wc will loosely refer to as fidelity. 



F 



„t(9-9)\ 



This name is suggested by the fact that Q < F < 1, with 
the values and 1 for the completely random guessing 
and perfect estimation, respectively. The relation with 
the Holevo variance can be read off from Eq. (1): F ^ 
(1 -I- y)~^/^. The average in F (or V) is over all possible 
outcomes as well as all possible signal states. 



F 



d9 



El uu 
I 9^ 



Me-e-^) 



27r 



tr[OxPe\ 



(3) 



where we have assumed a flat prior distribution on the 
circle. This is the natural choice when nothing is known 
beforehand about the phase we wish to estimate. An 
additional feature of the Holevo phase variance V is that 
it approaches the statistical variance (A^)^ = (6*^) — (6*)^ 
in the limit of accurate estimation (peaked distributions). 



i.e. 



V ~2{l-F)~ {Mf 



(4) 



To maximize F , it is useful to write the coherent state 
in the photon number eigenbasis (also referred to as Fock 
basis): pe = Y.ki^'^'"'''''^ Pki\k){l\- It is not difficult to 
prove that F is bounded by (see [22] for details) 



F = 



Y^e-^^^Y.P^^+^\-Ox\k+u 



<^|p,,+i|^|[0;,],+i,|, (5) 

k X 

which becomes an equality if [0^]fe+n: = |[0^]fe+n:|e'^'' 
for all k and x- Applying the Cauchy- Schwartz inequality 
and using positivity and the completeness of O^, which 
together amount to ^ |[0^]in,n| ^ 1, we finally obtain 
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(6) 



One can easily convince oneself that Holcvo's canonical 
phase measurement [15] 



[om, 



^-t4>(k-l) 



27r 



(7) 



saturates Eq. (6) and satisfies the unbiasedness condi- 
tion (2) [22]. Here, the matrix elements are again written 
in the Fock basis {|fc)}, </> is a uniform continuous param- 
eter (p e [0, 27r) that labels the outcomes, i.e., plays the 
role of X, and the optimal estimator is simply given by 
9^ = (/), with (/>£ [0,27r). 

At this point, the only task left is to compute the 
sum of matrix elements in Eq. (6). Recall that a co- 
herent state with mean photon number n is given by 



|a) = exp(-|a|V2)X;fca /Vfc!|fc), with |a| = ^. For 
a single copy of a coherent state it is straightforward to 
obtain from Eq. (6) that the maximum value of F is 
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It proves useful to cast the above expression into an inte- 
gral form, which is easier to study analytically. For this 
purpose, we use the identity 
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Vk+T 



dt 



-t{k+l) 



(9) 



and perform the (now trivial) summation over k in (8). 
By an appropriate change of variables one gets 



F 



dx- 



-xlal-^ 



v/-log(l-a;) 



(10) 



Eq. (10) is an integral representation of Eq. (8) and can 
be easily computed to arbitrary precision for any value 
of a. 

The large n regime can now be worked out. Prior to in- 
tegration, we Taylor-expand the logarithm around a; = 
in Eq. (10) and then integrate term by term. Note that 
each power of x that we retain gives a contribution of 
order l/|Qfp = 1/n. The upper limit of the correspond- 
ing integrals can be safely taken to be infinity, since this 
change will only contribute to differences that fall off ex- 
ponentially with n. To subleading order we obtain 



F = 1 



1 

8^ 



(11) 



This is the maximum average fidelity attained with a 
single, highly energetic coherent state {n ^ 1). Using 
Eq. (4) one has 
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(12) 



which agrees with the well known statistical variance of 
coherent states [23]. The alternative derivation we have 
presented here will prove very useful for SVS, as will be- 
come apparent in the remainder of the paper. 

Interestingly, the same analysis can be carried out for 
an arbitrary number of coherent states. The case of two 
copies N = 2 already contains all the ingredients of the 
solution for arbitrary N. We recall that a symmetric 
state with total photon number k is 
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simply realized, e.g., by a 50/50 beam splitter). To sim- 
plify the notation, we drop the mode labels {A, B, . . .) 
and the the vacuum states (|0)), as in (13) and (14), 
throughout this paper. Applying this reasoning induc- 
tively we see that an optimal generalized measurement 
on N identical coherent states \a)'^^ is formally equiva- 
lent to a single measurement on a coherent state \\/Na). 
Thus, for large N we get the HL relation 
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III. SQUEEZED VACUUM STATES 



We next address the case of SVS. Although the cal- 
culations are substantially more involved, the techniques 
are not so different from those presented in the preceding 
section. Here \i;g) is a SVS given by \ije) = U{e)S{r)\0), 
where S{r) is the squeezing operator S{r) ~ exp[r(a^^ — 
a^)/2], and U{6) is the phase shift operator already de- 
fined at the beginning of Section II. The parameter r is 
usually referred to as the squeezing parameter. The SVS 
in the Fock basis read 



l^«) = (l-/3^)^/'E 
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\2k), (16) 



where one can readily see that the SVS are superposi- 
tions of Fock states, |2fc), with an even number of pho- 
tons. In Eq. (16) we have defined /3 = tanhr, which 
in turn is related to the mean photon number through 

P=^n/{n+l). 

The A^-copy state vector |$e) = IV-'e)'^^ can be written 
as 



|$e) = (l-/?2)^/4E(/3e2'^)'=/ifc|2fc), (17) 

(18) 
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where |2fc) are the A^-mode symmetric states with total 
photon number 2k. E.g., for N ^ 2 they read 
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By using this definition we see that |a)®^ is unitarily 
equivalent to a coherent state with amplitude \^a. More 
precisely. 
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Hence, two identical coherent states |Q;)yi|a)B can be 
transformed by a two-mode unitary transformation into 
a single coherent state |-\/2a)c|0)n (this unitary can be 



where in this case /i^ = 1. The generalization for ar- 
bitrary A^ is straightforward. Note that the explicit 
form of hk guarantees that the normalization condition 
(<I>e|$e) = 1 is fulfilled. Actually, 
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(20) 



can be viewed as the probability mass function of a neg- 
ative binomial distribution [27] with failure probability 
given by /3^. 



The invariancc I'&e+Tr) = \^e) imposes a minor modi- 
fication of the Holevo phase variance (1) for SVS, which 
now reads 
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^2i{e-e)x 



1. 



(21) 



The factor of two in the exponent takes care of this in- 
variance at the expense of having a phase in the range 
[0, tt). Accordingly, the fidehty reads 



F=\{e^'^'-''^)l 



(22) 



and the relation between V and the statistical variance 
(A^)^ is now 



V ~2{l-F) ~A{Aef 



(23) 



From Eq. (22). and assuming that 9 is uniformly dis- 
tributed in the interval [0, tt), one can easily obtain that 
the bound to the fidelity is formally equivalent to Eq. (6): 
F < ^j, |pfefc-|_i|, where here p is written in the ba- 
sis of symmetrized states \2k) defined in Eq. (19) as 
Pkk' = {2k\p\2k'). Likewise a measurement that satu- 
rates this bound is given by 
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(also written in the same basis), and the optimal esti- 
mator is 9^ = 4>, with (j) S [0, tt). Taking into account 
Eq. (17), the explicit form of the bound reads 



F<(l-/32)^/2^/32'=+i/.,/.fc+i 



(25) 



k=0 



Using the techniques shown in Appendix A, this expres- 
sion can be recast into an integral form which is much 
easier to study, both analytically and numerically. It 
reads 



N 



F=il/3(l-/?2)^/2 



du'i 



^/.(^logu) , p,, 



(1 - u/32)^/2+l 



where Iq{x) is the zero order modified Bessel func- 
tion [24]. Eq. (26) follows from Eq. (A7) upon changing 
variables from r to u = (1 — r)//?^. 

Eq. (26) enables us to compute the Holevo variance in 
a very efficient way for arbitrary values of N and n. Fig. 1 
shows a log-log plot of NV against ?7 for 1 < A^ < 9. We 
see that for low n all curves have the same slope, while 
for large n, as we anticipated in the introduction, the 
slopes increase up to A^ = 4 and then stabilize. We next 
proceed to calculate these two scalings analytically. 



10 



0.1 - 






0.001 



10" 



10" 




0.01 



FIG. 1: Log-log plot of the scaled phase variance NV defined 
in Eq. (21) for A^ copies of SVS, 1 < AT < 9 (top to bottom), 
as a function of the mean photon number n (solid lines). For 
large n the lines become steeper as A'' increases, in agreement 
with Eq. (27). The slopes stabilize for N > 5. For small n the 
slopes are independent of A'^, as Eqs. (28) and (36) show. The 
dotted line is the limiting curve (N — » oo) for both small and 
large n. It follows from the HL (36), as discussed in Sec. V. 



this regime the crucial issue is to know the rate at which 
perfect estimation is achieved. For the sake of readabil- 
ity the technical details of the calculation are collected 
in Appendix A. In this section we only summarize and 
comment the main results. From Eqs. (A13) and (23) 
one readily sees that the statistical variance at leading 
order in 1/n is 
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(Af <3) 

{N = 4) (27) 

(^>5), 



where the values oiS,N ai'e given in Table I. We recall that 
these statistical variances are obtained assuming that one 
can perform the most general (collective) measurement 
on the A^-copy state pf . We find it remarkable that 
the 7i~2 scaling is only achieved for N > 5 (see Fig. 1). 
Notice that for a single copy, the optimal scaling is n~^/^, 
as compared to n"^^ which one would naively expect from 
HL. This scaling is even worse than that attained with 
coherent states n~^ [see Eq. (6)]. The same n~^^^ scal- 
ing shows up in the n-photon two mode state discussed 
in [25]. Using our analysis we can also study the optimal 
splitting of energy among copies. If one has a fixed, but 



A. Large squeezing 

For very energetic SVS, n ^ 1 (/3 ^- 1), the phase can 
be estimated with arbitrary accuracy; lim^^i F = 1. In 
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TABLE I: Values of the coefficients ^jv in Eq. (27). 
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FIG. 2: Log- log plot of the rescaled phase variance E V ^& 
a function of the total available energy E = nN (in units of 
Tiuj). The thin solid lines correspond (from top to bottom) 
to Ai" = 5, 6, 7; the thick green line to A*' = 8. Dashed lines 
correspond to A'^ > 9 in increasing order from bottom to top 
on the right side of the plot. For large E, variances scale as 
E^^ (all lines have horizontal asymptotes) and N = 8 clearly 
provides the smallest variance. 



large, amount of energy E that can be divided among N 
copies, each with mean photon number n = E/N, it is 
straightforward to obtain from Eq. (27) that asymptot- 
ically the optimal choice is A^ = 8. This result is also 
clear from Fig. 2. 



B. Small squeezing 

Let us now briefly focus on the low energy regime. This 
regime may be relevant in practical situations where the 
amount of squeezing is bounded by technological limita- 
tions. For /3 <C 1 we have /3 ~ ^/n and we can easily find 
the leading behavior of F by keeping the first terms in 
Eq. (25). For nN ^ 1, a simple calculation gives 



F : 
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(28) 



where the approximate expression for V follows from (21) 
and (22). This shows that the scaling is independent 
of A'', as is clear from Fig. 1. Note also that the fidelity 
only depends on the total energy E = nN. 



IV. INDIVIDUAL MEASUREMENTS: LARGE N 



tion between measurers dealing with the various copies. 
The Fisher information allows to address this question. 
In single-parameter estimation the Fisher information 
provides an asymptotic bound (the Cramer-Rao bound) 
for the accuracy that can be attained by repeating the 
same measurement on each copy of the sample. Recall 
that the Fisher information is defined as [28] 
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where p{(j)\6) is the conditional probability of obtaining 
the outcome (j) upon measuring on the state \tp0) that 
carries a phase 9 . In our case the measurement is defined 
by Eq. (24) and the corresponding outcome probabilities 
are 
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The covariancc of the measurement (24) implies that the 
Fisher information is independent of 0, as is apparent 
from the form of the conditional probability (30). For 
highly energetic SYS (/? ~ 1) we get 
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(we refer to Appendix B for details), and by using the 
Cramer-Rao bound N{A9)'^ > I~^, we obtain 
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Recall that this bound is attainable for a large sample, 
N -^ oo, with, e.g., a Bayesian or maximum likelihood 
estimator [18]. 

Hence, we obtain that the optimal measurement per- 
formed sequentially on a large sample of identically pre- 
pared SVS gives an accuracy that scales as the HL, a 
result that could not easily be anticipated. Notice, how- 
ever, that the coefficient of the variance is more than a 
factor of two larger than that of the HL (7/24 as com- 
pared to 1/8). This means that in order to obtain the 
optimal accuracy we may require classical communica- 
tion [14], i.e., adaptivity of the measurements. 



We have shown that the accuracy that can be achieved 
with a single SVS increases with n at a, much lower rate 
than that inferred from the HL. The natural question 
that arises is whether the latter can be achieved by per- 
forming the same measurement on a large set of identical 
copies. In other words, we wish to know if the attainabil- 
ity of the HL requires some sort of classical communica- 



V. DISCUSSION 

From Eq. (27) and Fig. 1 we see that the HL, which pre- 
dicts a scaling n~^ for the variance, is only achieved for 
N > 5. Had we naively extrapolated the scaling n~^/^ 
for A^ < 4, we would have predicted a breakdown of the 
HL (see Appendix A). Of course, this is not the case: for 



A^ > 4 the terms with the dangerous exponent behavior 
in Eq. (A7) become sub-dominant in the large n hmit, 
and the statistical variance scales as dictated by the HL. 
We can also obtain the exact dependence in the mean 
number of photons for large N. The easiest way is to 
rewrite hkhk+i in Eq. (25) as hl^y{N/2 + k)/{k + 1), so 
that 
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where the negative binomial probability distribution pk 
is defined in Eq. (20). To compute this expectation value 
we expand the square root around the mean of the distri- 
bution, (k) = {N/2)P^/{1 - /32), up to second order. By 
recalling that (Afc)^ = (fc)/(l — /3^), it is straightforward 
to obtain 
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(35) 



or equivalently 



performed on each copy of the sample. It gives rise to 
a scaling of the variance which agrees with the HL of 
SVS up to a constant multiplicative factor of the order of 
two. In some sense, this result relaxes the need for adap- 
tive protocols at the expense of having an asymptotic 
rate that is roughly a half that of the optimal protocol. 
This shows that individual non-adaptive measurements 
can harness the enhanced phase variance offered by SVS. 



VI. ACKNOWLEDGMENTS 

We are grateful to A. Acin, J. Calsamiglia, M. Mitchell 
and M.B. Plcnio for very useful discussions, and to 
D. Diego, for contributing to the very early stages of 
this work. We acknowledge financial support from the 
Spanish MCyT, CONSOLIDER2006-00019 and FIS2005- 
01369, and CIRIT 2005SGR-00994. A.M. also acknowl- 
edges financial support from QIP IRC GR/S82176/01. 



APPENDIX A: CALCULATION OF F FOR SVS 
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(36) 



which is the exact expression of the HL (see [14, 26]^ and 
references therein). 

Our results have a number of implications both of 
fundamental and of practical interest. In the high en- 
ergy regime they show that, somewhat unexpectedly, the 
phase resolution of a single squeezed vacuum state is 
worse than that of a coherent state with the same energy. 
Our results indicate that whenever a phase measurement 
is to be performed through interaction with a single copy 
(or mode) of the probe, one is better off using coherent 
states. This scenario changes significantly as one moves 
to the multi-copy case. For two copies the resolution is 
already comparable to that of coherent states and the 
variance reaches the HL scaling (n~^) for A^ > 4. It is 
important to notice that we have considered arbitrary 
collective measurements and, hence, our results give the 
ultimate precision bounds allowed by quantum mechan- 
ics. In addition, our analysis provides the means to de- 
termine the optimal energy splitting among copies, show- 
ing that for highly energetic SVS, the optimal choice is 
A^ = 8, as shown in Fig. 2. 

In the low energy regime we have shown that the vari- 
ance (fidelity) is only a function of the total available 
energy, regardless the way it splits among copies. This is 
relevant for practical implementations of squeezed state 
metrology, where usually the amount of available squeez- 
ing is small. 

We have also analyzed the asymptotic accuracy of 
Holevo's canonical phase measurement (24) when it is 



1 Note the wrong sign in Eq. (10.39) of [26] 



In this Appendix we provide details of our calculation 
of the fidelity for SVS. Our starting point is Eq. (25). 
Taking into account Eq. (18) we can write 



F<{l-(3^f^^J2 



a' 



•2k+l 



Ck 



-i^{k + l)ik + N/2) 



fe=0 



where we have defined 



Cfc 



N fk + N/2 



(Al) 



(A2) 



We next apply the integral representation (9) to both 

(fc + l)-i/2 and {k + N/2)-^^^, 



1 



^{k + l){k + N/2) 
Ja TTy/trh 



(A3) 



This enables us to sum up the series (over fc) in (Al) 
by recalling the negative binomial expansion ^^ c^z ~ 
{N/2){1 - z)-^/2-i. We are left with a double integral, 
which can be further simplified by the change of vari- 
ables (ti,t2) -^ (u,t): 



t2 = -(l-w)log 



1-T 



(A4) 

(A5) 



and by recalling the integral representation of the zero 
order modified Bessel function of the first kind [24] 



1 [^ f,{l-2u)x 

Io{x) = — / du 



TT Jo V"(l-") 



(A6) 



The remaining integral can be cast as 
where 



(A7) 



where ^at arc formally defined in terms of definite inte- 
grals. Its numerical values are given in Table I. The exact 
value ^2 = 1/2 can be trivially obtained from Eq. (34), 
which yields F = j3. 



9ir,P) 



1 (\-T 



\-t\ 13' 



/o(2^1og^ 



(A8) 



Eqs. (A7) and (A8) provide a very useful expression 
of F. The integral over r can be computed to arbitrary 
accuracy and it is valid for any number of copies and for 
any average photon number. It allows, e.g., to find out 
the optimal splitting of the available energy E among 
copies for phase estimation (see Sec. Ill A and Fig. 2). 

We next derive from ( A7) the high energy scaling of the 
fidelity. Let us first introduce some short-hand notation. 
By N{P), we denote the prefactor in (A7), 



TV 



AA(/3) = -/3(1 - /32)^/2, 



(A9) 



and notice that it is of order n ^/^. We further define 
Sp{t,P) = Y.k=o^k{l3)T'^ to be the truncated Taylor 
expansion of g{T,(3) around r = up to order P. One 
can easily check that its corresponding integral over r. 



Qm 



S[N/2] {t, P) 



1-/3^ 



-1 + 7V/2 



dr. 



(AlO) 



is of order n^''^ {[x] is the integer part oi x), whereas the 
integral of the remainder, 



52m = 



girJJ) -S'[Ar/2](r,/3) 



1-/32 



-l+N/2 



dr. 



(All) 



is of order n'^ . The leading contribution can be computed 
by taking the limit (3 -^ I in Q2- Thus, up to subleading 
order, we have 



F = Af{p)[gm + 52(1)]. 



(A12) 



We next expand the first term, J\f{(3)Qi{f]), in (A12) in 
powers of n~^'^. This expansion provides the leading 
order contribution to the fidelity (the unity), a sublead- 
ing term of order n^'^, which becomes dominant only for 
A'^ > 4, and a series of inverse power contributions start- 
ing at order 7i~^'^. These latter become irrelevant if 
A^ > 4. On the other hand, the second term, J\f{/3)G2{i), 
results in corrections of order n~ " which contribute 
to the subleading order if A^ < 4 and become irrelevant 
otherwise (of order smaller than n~^). 
In summary, we have 



r e 



N 



F^l- ■( 



nN/2 
(l/8)logn + e4 

1 

_ 4(A - 4)n2 



{N<3) 

{N = 4) (A13) 

(A^ > 5), 



APPENDIX B: PHASE ESTIMATION WITH 
INDIVIDUAL MEASUREMENTS 

In this Appendix we compute the Fisher information 
of the optimal 1-copy measurement. This leads to the 
Cramer-Rao bound discussed in section IV. 

We first notice that the dominant behavior of S{6) in 
Eq. (31) is determined by the asymptotic expressions of 
the factorials in the sum. The Stirling approximation 
gives 



SW 



P' 



1-13' 



1/4 



1/4 Z^ 



2i(p\k 



/ci/4 



1/4 



Lii/4(/3i 



,2i4,\ 



ri/4 



(Bl) 



where Lis(z) = J^kLi ^'^ / ^"^ i^ ^^^ polylogarithm func- 
tion of order s [24] . Since we are only interested in the 
asymptotic behavior for large squeezing {f3 ^ 1), we can 
use the first order expansion around z = 1, which reads 
Lii/4(z) = r(3/4)(l - z)-3/4 + . . . [24]. We note in pass- 
ing that the probability law obtained by retaining only 
this term is consistent with the normalization condition, 
i.e., J p{(j)\9)d4> — > 1 when /3 ^ 1. In this limit we can 
easily compute the derivatives required in (29) and ob- 
tain 



9r(|)^yr^^ 

7r3/2 



sin^ 20 



(l-2/3cos20-H/32)- 



(B2) 



where we have used the rotational invariance of the in- 
tegral. Notice that the Fisher information is indepen- 
dent of the phase of \tpg). This is just a consequence 
of the continuous and covariant character of the mea- 
surement (7) and the isotropy of the prior. Again, in 
the limit (3^1 the main contribution to the integral 
is peaked around 0^0 and we can Taylor-expand the 
trigonometric functions and safely extend the integration 
limits from (— 7r/2,7r/2) to (— oo,cx)). We thus obtain 



36r^(|)/3^v/r^ /■-_ 

^3/2 loo [(1 



3 V2/3(l + /3) 
7 (1-/3)2 



-00 [(l-/3)2+4/?(/.2]T 

6 



7(1-/3)2 



(B3) 



This is the Fisher information used in Eq. (32) in the 
main text. 
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